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MEASUREMENT OF AERODYNAMIC FORCES FOR VARIOUS MEAN ANGLES OF ATTACK ON AN 
AIRFOIL OSCILLATING IN PITCH AND ON TWO FINITE-SPAN WINGS OSCILLATING IN 
BENDING WITH EMPHASIS ON DAMPING IN THE STALL 1 

By A. Gerald Rainey 


SUMMARY 

The oscillating air forces on a two-dimensional wing oscil- 
lating in pitch about the mid chord have been measured at various 
mean angles of attack and at Mach numbers of 0.35 and 0.7. 
The magnitudes of normal-force and pitching-moment coefficients 
were much higher at high angles of attack than at low angles of 
attack for some conditions. Large regions of negative damping 
in pitch were found, and it was shown that the effect of increasing 
the Mach number from 0.35 to 0.7 was to decrease the initial 
angle of attack at which negative damping occurred. 

Measurements of the aerodynamic damping of a 10-percent- 
thick and of a 3 -percent-thick finite-span wing oscillating in 
the first bending mode indicate no regions o f negative damping 
for this type o f motion over the range of variables covered. The 
damping measured at high angles o f attack was generally larger 
than that at low angles of attack. 

INTRODUCTION 

Recently, the problems of stall flutter and buffeting have 
been the subjects of several analytical investigations. (See 
refs. 1 to 8.) These studies have indicated some interesting 
approaches to the dynamic problems associated with sep- 
arated flow; however, because of the great difficulty in 
treating the unsteady aerodynamics of wings with separated 
flow, a completely satisfactory analytical approach has not 
been developed. Therefore, additional experimental infor- 
mation is required in order to describe further the phenomena 
of stall flutter and buffeting. 

It is generally accepted that, when a wing is sufficiently 
stalled, fluctuating forces act on the wing and, if the wing is 
flexible, it will respond to these forces. If the total damping 
of the system remains positive, the characteristics of the 
response will depend in part on the characteristics of the 
fluctuating forcing function, and the amount the wing 
responds will depend on the relative magnitudes of the 
forcing function and the damping. This phenomenon may 
serve as a rather general description of buffeting. 

If, on the other hand, the damping becomes negative, the 
characteristics and amount of the response will depend 
primarily on the character of the damping, that is, on the 
degree of nonlinearity. This latter case describes the phe- 
nomenon of stall flutter. Then, in essence, the phenomena 


of stall 11 ut ter and buffeting, as defined herein, differ primarily 
in the sign of the damping. Consequently, knowledge of 
the conditions under which the aerodynamic damping 
becomes negative is essential in order to decide on the proper 
approach to the problem. 

The purpose of this report is to present the results of two 
somewhat different investigations both of which were in- 
tended to provide additional information concerning the 
conditions under which aerodynamic damping may become 
negative, as well as some quantitative information concerning 
the amount of aerodynamic damping. In addition, the 
techniques and instrumentation employed are given in some 
detail. The first part of this report deals with the results 
of some measurements of oscillating normal forces and 
pitching moments for an NACA 65 AO 10 airfoil oscillating 
in pitch about the midchord. Some additional results 
obtained with this pitching airfoil with two different types 
of disturbances at the leading edge are presented in appendix 
A. The second part gives the results of measurements of 
aerodynamic damping for two different finite-span wings 
oscillating in essentially the first bending mode. 

SYMBOLS 

a x tip deflection in first bending mode, ft 

c wing chord, ft 

c T reference wing chord at three-quarter semispan, ft 

C structural damping per unit length, lb-sec/ft 2 

C cr critical value of damping, lb-sec/ft 

C a effective value of aerodynamic damping, lb-sec/ft 

C s effective value of structural damping, lb-sec/ft 

c m section pitching-moment coefficient about mid- 

chord 

dc m absolute magnitude of fundamental component of 

da oscillating section pitching-moment coefficient 

about midchord per radian 
c n section normal-force coefficient 

dc n absolute magnitude of fundamental component of 

da oscillating section normal-force coefficient per 

radian 

El bending rigidity, lb-ft 2 

/ frequency, cps 

F total applied force acting on finite-span wing, lb 


1 Supersedes NACA Technical Note 3643 by A. Gerald Rainey, 1956. 
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effective value of aerodynamic-damping coefficient 
for finite-span wing oscillating in bending mode 
(eq. (3)) 

vertical deflection of equivalent spring-mass sys- 
tem., ft 

first-bending-mode shape expressed in terms of 
unit tip deflection 
reduced frequency, uc/2V 

effective spring constant of wire, coil, and spring 
attached to finite-span models, lb/ft 
spring constant of equivalent spring-mass system, 
lb/ft 

beam stiffness in first bending mode, lb/ft 
length of semispan models, ft 
mass per unit length, lb-sec 2 /ft 2 
absolute magnitude of fundamental component of 
oscillatory aerodynamic pitching moment about 
midchord, ft-lb 

effective mass of wire, eoil, and spring attached to 
finite-span models, lb-sec 2 /ft 
increment of coil mass, lb-sec 2 /ft 
mass of equivalent spring-mass system, lb-sec 2 /ft 
beam mass in first bending mode, lb-sec 2 /ft 
total oscillatory pitching moment indicated by 
pressure-gage integrator, ft-lb 
part of pitching-moment signal passed by filter 
which is definitely related in phase to the 
motion, ft-lb 

part of pitching-moment signal passed by filter 
which is not correlated with the motion, ft-lb 
aerodynamic and exciting forces normal to finite- 
span wings 

absolute magnitude of fundamental component of 
oscillatory aerodynamic normal force, lb 
total oscillatory normal force indicated by pres- 
sure-gage integrator, lb 

part of normal-force signal passed by filter which 
is directly related in phase to the motion, lb 
part of normal-force signal passed by filter which 
is correlated with the motion, lb 
incremental pressure acting between upper and 
lower surface of wing, lb/sq ft 
wing area, sq ft 
arbitrary period of time, sec 
time, sec 

flow velocity, ft/sec 
spanwise coordinate, ft 

spamvise coordinate at point of attachment of 
wire-spring-shaker-coil system, ft 
instantaneous angle of attack, radians 
absolute magnitude of fundamental component 
of oscillatory angle of attack, deg (except when 
used in defining oscillatory derivatives) 
mean angle of attack about which oscillation 
occurs, deg 


q?oi part of angle-of-attaek signal passed by filter 

which is directly related in phase to normal 
force or pi telling moment 
K e effective value of mass-ratio parameter 

p mass density of test medium, lb-sec 2 /ft 4 

0 m phase angle by which pitching moment leads the 

motion, deg 

cj) n phase angle by which normal force leads the 

motion, deg 

co circular frequency, 2wf , radians/sec 

co! first natural frequency, radians/sec 

co lA first natural frequency of wing-shaker system with 

added mass, radians/sec 

Subscripts: 

0 quantities which have same frequency and which 

are directly related in phase t o sinusoidal forces 
or motions 

R quantities which are not correlated with sinusoidal 

forces or motions 

Bars indicate an average with respect to time. 

One dot indicates first derivative with respect to time, 
and two dots indicate second derivative with respect to 
time. 

TWO-DIMENSIONAL AIRFOIL OSCILLATING IN 
PITCH ABOUT MIDCHORD 

TECHNIQUE 

The normal forces and pitching moments acting on an 
NACA 65A010 airfoil oscillating in pitch about the midchord 
at high mean angles of attack were measured by use of a 
technique involving an electrical pressure-gage integrator. 
Some of the methods and instrumentation used have been 
briefly described in reference 9 wherein results are presented 
of measurements of oscillating air forces at a mean angle 
of attack of zero. For this condition, the instrumentation 
and methods presented in reference 9 were considered to 
be adequate; but, when measurements were undertaken at 
mean angles of attack beyond the stall, it was found that 
relatively large random forces existed and somewhat different 
instrumentation and data-reduction procedures were required 
in order to obtain the desired quantities. A detailed de- 
scription of the apparatus, testing procedures, and methods 
of reducing the data follows. 

DESCRIPTION OF APPARATUS 

Wind tunnel. — The wind tunnel used in this investigation 
was the Langley 2- by 4-foot flutter research tunnel which 
is a single-return type equipped to use either air or Freon-12 
as a test medium at pressures from 1 atmosphere down to 
about Ys atmosphere. The actual pressure and medium used 
will be discussed subsequently. 

Airfoil and oscillator. — The airfoil used in determining the 
oscillatory normal forces and moments due to pitch had an 
NACA 65A010 airfoil section and a constant chord of 1 foot. 
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In addition to the tests of the clean airfoil, a few limited 
tests were performed with two different types of leading- 
edge disturbances and these are presented in appendix A. 
The model completely spanned the 2-foot dimension of the 
test section and was equipped with end plates which rotated 
with the airfoil to prevent tip leakage. The model was 
constructed from aluminum alloy that was machined in the 
form of two half shells with a solid section in the center to which 
were attached the pressure gages and the position indicator. 
The construction of the model is illustrated in figure 1 . 

The airfoil oscillator and general test configuration is 
illustrated in figure 2. Power to drive the wing is furnished 
by a synchronous-drive motor with a variable-frequency 
voltage supply. The speed of the drive motor, which was 
continuously controllable from 300 to 3,600 rpm, determined 
the frequency of oscillation. Attached to the motor shaft 
on both sides of the test section is the inner race of a bearing 
called an eccentric cam in figure 2. The eccentricity of this 
cam determines the amplitude of the oscillation. The 
oscillation is transmitted to the wing through a drive column 



Figure 2. — Schematic drawing of model and oscillating mechanism. 



Strain gages 

Torque rod 


Pressure - gage cavities 


Air cavity 


Pressure gage 


Orifice | 
locations 



Percent 2.5 5.0 10.0 15.0 20.0 

chord 


Orifice 
30.0 40.0 


-Neoprene seal 

50.0 60.0 


70.0 


80.0 


90.0 

L- 91763 

Figure 1 . — Exploded view and cross section of model showing location of pressure gages, orifices, and position indicator. 
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(equipped with a turnbuckle for adjustment of mean angle 
of attack) pinned to a rocker arm which rotates a drive shaft 
attached to the airfoil. 

Instrumentation. — The differential pressure acting on the 
wing was converted to an electrical signal by use of NACA 
miniature electrical pressure gages model 49-TP. (See ref. 
10.) The 12 ’gages were located at 2.5, 5, 10, 15, 20, 30, 40, 
50, 60, 70, 80, and 90 percent chord. The pressure-gage 
installation is illustrated in figure 1. The gages between 5 
and 80 percent chord were clamped between the two wing 
halves and were connected to 0.030-inch-diameter orifices 
by 0.0625-inch-diameter drilled holes which had depths that 
varied with the thickness of the wing. The gages at 2.5 and 
90 percent chord were connected to orifices by tubing which 
had an inside diameter of 0.040 inch and had lengths of 
0.6 and 1.3 inches, respectively. The experimentally deter- 
mined dynamic-response curves for the gages at 2.5, 90, and 
20 percent chord are shown in figure 3. The resonant peaks 
occurring at high frequencies are associated with an acoustic 
phenomenon. The response curves presented in figure 3 
were obtained at atmospheric pressure in air. The response 
characteristics of the pressure-gage-orifice installations have 
not been determined with Freon-12 as the ambient medium; 


however, it is believed that the systems have fiat response 
to a sufficiently high frequency so that no response errors 
were involved in the measurements contained in this investi- 
gation which covers the frequency range from 0 to 35 cps. 

The pressure-gage outputs were weighted and integrated 
electrically to give a signal proportional to lift or pitching 
moment. The weighting was computed on the basis that 
each gage represented the area or distance halfway between 
adjacent gages. For example, the first gage (at 2.5 percent c) 
was assumed to represent the area from the leading edge to 
the point halfway between the first gage and the second gage 
or an area of 3.75 percent of the total area. The second 
gage was assumed to represent the area or distance between 
3.75 percent c and 7.5 percent c or an area of 3.75 percent 
of the total area. For pitching moments, each area was 
multiplied by the distance from its centroid to the refer- 
ence axis. The principles and details of electrical integration 
are discussed in reference 1 1 . 

The angular motion of the center-line station of the airfoil 
was converted to an electrical signal by means of a torsion 
rod which is illustrated in figure 1. A resistance-wire- 
strain-gage bridge arranged to indicate torsional strains was 
attached to a steel rod 0.5 inch in diameter which had one 



Figure 3.— Dynamic response characteristics of three pressure-gage-orifice installations. 
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end fixed to the center of the airfoil and the other end fixed 
to the tunnel wall. The first natural frequency of the rod 
in torsion was estimated to be about 700 cps. 

A block diagram of the instrumentation is shown in 
figure 4. The electrical signals proportional to normal 
force or pitching moment and angular displacement were 
operated on in a variety of ways by suitable switching 
arrangements. 

Time-history records of the lift or pitching moment and 
angular displacement were obtained by use of a recording 
oscillograph. Some observations on the time histories of 
individual pressures on an oscillating airfoil near the stall 
are presented in appendix B. 

The magnitude of the fundamental normal -force or pitch- 
ing-moment vector was determined by feeding the signal 
through a narrow band-pass variable-frequency filter which 
was tuned to the frequency of oscillation. The signal out 
of the filter was then amplified through a linear amplifier 
and fed into a vacuum thermocouple, the output of which 
was read on a heavily damped microammeter. This heavily 
damped mean-square indicating device was required because 
of the fluctuations in normal force or pitching moment 
encountered at the higli-angle-of-attack conditions. 

The phase angle of the fundamental component of the 
normal force or pitching moment with respect to the angular 
displacement was obtained bv passing both the normal-force 
or moment and angular-displacement signals through identi- 



4>n or <t>m and " 

Figure 4. — Block diagram of instrumentation used for two- 
dimensional tests. 


cal channels of the narrow band-pass filter then into identical 
channels of a pulse-shaping device. The pulse shaper con- 
sisted, essentially, of several stages of amplification, clipping, 
and differentiation which resulted in conversion of the input 
sine wave to a series of sharp pulses corresponding in time 
to the positive crossover points of the sine wave. The 
phase angle between normal force and position was then 
obtained by starting an electronic chronograph with the 
leading signal and stopping the chronograph with the lagging 
signal. The time obtained in this manner was the time lag 
between normal force or moment and position which, when 
divided by the period of the oscillation and multiplied by 
360°, yielded the phase angle in degrees. A similar method 
of measuring phase angles, which was developed independ- 
ently, is described in detail in reference 12. The period of 
the oscillation was obtained by starting and stopping the 
chronograph with a single signal, usually the angular dis- 
placement. Of course, slight differences in components and 
tuning caused extraneous shifts in phase between the two 
signals. This tare phase angle was measured during each 
measurement of phase angle by feeding a single signal 
(usually the angular displacement) through both channels. 

PROCEDURE 

The aerodynamic normal forces and pitching moments 
acting on the two-dimensional airfoil oscillating in pitch were 
determined at Mach numbers of 0.35 and 0.7 and at a 
Reynolds number of 5.3X10 6 . Freon-12 was used as a test 
medium. The Reynolds number was held constant by 
reducing the tunnel pressure from 1 atmosphere for the tests 
at a Mach number of 0.35 to about V 2 atmosphere for the 
tests at a Mach number of 0.7. 

After the desired tunnel velocity and pressure had been 
obtained, the mean angle of attack was adjusted by use of 
the turnbuckles shown in figure 2. The zero-frequency 
increment in normal force or moment was then obtained with 
the model set at its two extremes; this setting was accom- 
plished by turning the oscillator shaft through one complete 
revolution very slowly. The incremental voltage corre- 
sponding to the incremental force or moment was obtained 
from a self-balancing potentiometer. 

After the zero-frequency measurements were made, the 
drive motor was started and brought to a speed correspond- 
ing to an airfoil oscillation frequency of approximately 10 
cps. The lower limit of frequency was determined by the 
low-frequenc.y limit of the narrow band-pass variable- 
frequency filter. After the desired frequency had been 
reached, the normal-force or pitching-moment and position 
signals were operated on in the manner previously described 
to obtain the quantities of interest. The frequency was 
then increased in increments of about 5 cps to a maximum 
frequency of about 35 cps. Data were not obtained at 
frequencies above 35 cps because of the uncertainties of the 
effect of wind-tunnel walls as indicated in reference 9. 
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This procedure was repeated at initial angles of attack 
of 0°, 4°, 8°, 12°, and 16° at a Mach number of 0.35; also 
0°, 4°, 6°, and 8° at a Mach number of 0.7. The amplitude 
of oscillation for all conditions was about 1.2°. 

DATA REDUCTION 

The forces acting on an oscillating wing under separated- 
flow conditions may be thought of as being composed of 
three parts: (1) the steady or mean component, (2) those 
forces which are related in time to the motion, that is, have 
a definite phase relationship, and (3) those forces which are 
not correlated with the motion. The third type of unsteady 
forces will be random in character and probably will be 
similar, at least, to the fluctuating forces which act on a 
stationary rigid wing when it is stalled. At low angles of 
attack, a random force existed which was small (about twice 
the instrument noise level) and, presumably, was due to 
residual tunnel turbulence. At high angles of attack, the 
random forces were of the order of ten times the random 
forces which existed at low angles of attack; and, conse- 
quently, consideration had to be given to ways and means of 
accounting for these forces in making the desired measure- 
ments. 

This composition of forces may be expressed simply as 
follows: 

N(t)=N+N 0 (t)+N R (t) 

where N denotes the steady-state force associated with the 
mean angle of attack a, N 0 (t) refers to the oscillating forces 
which are related to the sinusoidal motion a(t)=\a\e iut , and 
N R (t ) designates the random forces which are not correlated 
with the motion. The problem, then, is to separate properly 
the three types of forces so that the significant properties 
can be determined. This separation can be done in several 
ways. The method used in these tests of the two-dimen- 
sional model consisted of passing the normal-force or pitch- 
ing-moment signal through a narrow band-pass variable- 
frequency filter tuned to the frequency of oscillation. In 
this manner, a large part of the random forces, as well as 
the mean force, were rejected; therefore, only the fundamen- 
tal force or moment vector and the random parts lying 
within the band pass of the filter were left. The possible 
existence and importance of higher harmonics in the force or 
moment is discussed in appendix C. 

The output of the filter is illustrated vectoriallv in figure 5. 
The large vector representing the fundamental component 
of normal force or pitching moment is fixed in phase, but 
the small vector representing the random part of the force 
or moment not rejected by the filter rotates randomly about 
the end of the fundamental component and varies randomly 
in length. The resultant vector will vary randomly in both 
magnitude and phase over a small range. Since the phase- 
measuring technique yielded the phase angle for a single 
cycle of oscillation, it w^as necessary to obtain several values 
for a given condition. The arithmetic mean of the several 



Figure 5. — Vector diagram of filter output tuned to frequency of 

oscillation. 


values was used as the best value to represent the phase of 
the fundamental component. The number of values taken 
at each test condition was increased or decreased according 
to the range of fluctuation encountered. Usually about 10 
values were considered sufficient. 

The magnitude of the fundamental component of normal 
force or pitching moment was obtained by passing the output 
of the filter through a mean-square indicating device which 
performed the following operation: 

[v 01 (o+A^(Ol 2 ^=V^+2iV^+Aiv7 2 

where iVoi is the fundamental component, A N R is the part 
of the random force lying within the band pass of the filter, 
and T is an arbitrary period of time. Since the time average 
of the product of two uncorrelated quantities must be zero, 
the reading of the mean-square meter will be proportional to 
N 01 2 +AN r 2 . The value of AY was obtained by subtracting 
the meter reading with the airfoil held rigid at the mean 
angle of attack from the meter reading during the oscillation. 

This operation assumes that the random forces are not 
affected by the oscillation. The severity of this assumption 
involves at least two conditions. Since the intensity of the 
random forces is a function of angle of attack, a modulation 
of the random forces at the frequency of the variation 
in angle of attack might be expected. This effect has been 
minimized by obtaining the estimates of AN r 2 at the mean 
angle of attack. 

An additional effect of the oscillation on the random forces 
might result from a change in the general character of the 
separated flow from the stationary case to the oscillating 
case. For example, the oscillation of the airfoil might alter 
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the boundary-layer stability in such a way as to cause a 
change in both the intensity and frequency content of the 
fluctuating forces. The possibility of a flow change of this 
type is a difficult one to check; however, some information 
supporting the validity of the assumption used in obtaining 
N 0 { 2 from the meter readings has been obtained. With the 
airfoil held stationary at the mean angle of attack; an esti- 
mate of the frequency content of the random forces was 
obtained by tuning the narrow band-pass filter to various 
frequencies within the range of interest. Similar measure- 
ments were made with the airfoil oscillating at a particular 
frequency about the mean angle of attack. The frequency 
content for these two conditions are compared witli the 
filter characteristic function in figure 6. The data for the 
condition where the wing was held stationary indicate a 
relatively smooth, low level of fluctuating forces. The 


data obtained with the wing oscillating are similar to that 
for the stationary wing except in the vicinity of the oscillation 
frequency. Near the wing oscillation frequency, the curve 
representing the filter characteristic function would be a good 
fairing of the data obtained with the wing oscillating; and, 
since the filter characteristic curve represents the filter 
output to a single frequency input, it can be concluded that 
the primary effect of oscillating the airfoil on the frequency 
content of the air forces was the addition of a single com- 
ponent at the frequency of oscillation. 

In addition to the aerodynamic normal force, the vector 
No contains a small component due to the inertia effects 
on the pressure gages. The magnitude of this tare-accelera- 
tion effect was determined by oscillating the airfoil in a near 
vacuum and is shown in figure 7. These tare values were 
subtracted vectorially from the N ol determinations. The 



Figure 6. — Comparison of normal-force filter output with and without oscillation with filter characteristic function. a=12° ; /=35 cps; 

and Mach number of 0.35. 
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Figure 7. — Apparent normal force and pitching moment due to inertia of pressure-gage diaphragm. 


tare values were kept small relative to the aerodynamic 
components by conducting the tests at as high a value of 
dynamic pressure as was consistent with the required Mach 
number and Reynolds number. Lines representing 1 percent 
of the smallest oscillating normal force and pitching moment 
are indicated in figure 7. The tare value in the worst case 
amounted to about 4 percent of the magnitude of the 
fundamental component of pitching moment and caused a 
phase correction of about 2°. 

Because of the inertia and flexibility of the mechanical 
system, the absolute magnitude of the pitching oscillation 
was a function of frequency even though the motion of the 
positive-displacement drive mechanism remained at the 
design value of |a| = 1.2°. It was found that the airfoil 
itself did not distort appreciably along its span, but there 
w r as a sizeable twisting of the drive shafts relative to the 
rocker arms outside the test section. Tests indicated that 


this twisting or distortion of the system was not affected 
by air forces; consequently, it w r as possible to use a single 
function to represent the variation of the absolute magnitude 
of the pitching oscillation with frequency. This function is 
shown in figure 8. 

After the absolute magnitudes of the fundamental com- 
ponents of aerodynamic normal force and pitching moment 
had been determined from the meter readings by use of 
suitable calibration constants, they were converted into 
coefficient form according to the following definitions: 


dc„ 

|A T | 

da 

M 

dc m 

\M\ 

da 

lf>V>Sc\a 
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Figure 8. — Variation with frequency of absolute amplitude of 
oscillation of two-dimensional airfoil. 


where |a| is measured in radians. 

The sign convention used in orienting the vector defines 
upward normal force, nose-up pitching moment, and nose-up 
pitch deflection as positive. 

The component of pitching-moment coefficient which is in 
phase with the angular velocity is related to the damping; 

dc 

therefore, the damping coefficient is defined as n 


da 


sm </>„ 


The sign convention used in describing the coefficients of 
damping in pitch is illustrated in the following vector 
diagram: 


Positive or stable damping Negative or unstable damping 




dc m 

da. V/ 
/9m 

dc m 

da Sm 


1 da \ 

dc m 



da 


RESULTS AND DISCUSSION 


The magnitude and phase of the fundamental components 
of oscillatory normal-force and pitching-moment coefficients 
for the two-dimensional NACA 65 AO 10 airfoil oscillating in 
pitch about the midchord are shown in figures 9 to 12 as 
functions of reduced frequency for various angles of attack. 


The coefficients of damping in pitch 


dc„ 


sin 4> m are shown as 


functions of reduced frequency in figures 13 and 14. 


LOW ANGLES OF ATTACK 

For purposes of comparison, the theoretical normal forces 
and pitching moments calculated for compressible, two- 


dimensional flow are shown along with the measured two- 
dimensional data in figures 9 to 12. The theoretical coeffi- 
cients used are those which are presented in reference 13 and 
whfcli are modified to agree with the sign convention and 
coefficient form used in this report. 

Examination of figures 9 and 10 indicates that the forces 
and moments measured at a = 0° and a=4° compare very 
well with the theoretical values at a Mach number of 0.35. 
The conventional steady-state section normal-force and 
section pitching-moment coefficients have been measured for 
a Mach number of 0.35 and a Reynolds number of 5.3 X10 6 
and are presented as functions of angle of attack in figure 15. 
These data indicate a linear relation between both normal- 
force and pitching-moment coefficients and angle of attack 
up to an angle of attack of about 7°. Tuft studies conducted 
at these same conditions indicated no disturbances on the 
upper surface until an angle of attack of 6° was reached. 
With the airfoil in a clean condition (that is, with no tufts 
attached), oscillographic observation of the outputs of the 
individual pressure gages indicated no pressure fluctuations 
up to an angle of attack of 7.6°. It appears reasonable, 
then, that there were no important boundary-layer effects on 
the static normal forces and pitching moments for angles of 
attack up to at least 6°. Furthermore, the oscillatory 
derivatives are not expected to vary appreciably with the 
mean angle of attack so long as the airfoil remains in the 
unseparated-flow condition throughout the range of oscilla- 
tion . 

Examination of figures 11 and 12 indicates that, at a 
Mach number of 0.7, the data for both a = 0° and a = 4° differ 
appreciably from theory. The discrepancy between theory 
and experiment at <x = 0° may be attributed primarily to 
effect of wind-tunnel walls. No corrections have been 
applied to the data to account for the presence of the walls. 
Botli the theory and experiments of reference 9 indicate an 
increase in the effect of tunnel walls when the Mach number 
is increased from 0.35 to 0.7. (Ref. 9 gives a comparison 
between theory and experiments on the same two-dimen- 
sional model as that used in the present investigation. This 
comparison shows good agreement between the measured 
forces and moments and those calculated which included the 
effect of the walls.) 

The discrepancies between the data for a =4° and the 
theory may be due, in part, to tunnel-wall effects or to 
boundary-layer and shock effects. There is virtually no 
information available concerning the effect of angle of attack 
on the wind-tunnel-wall corrections for oscillating air forces ; 
however, it appears reasonable that the effects of the walls 
might be greater at higher angles of attack because of the 
increased blockage. The static normal-force and pitching- 
moment coefficients measured at a Mach number of 0.7 
are shown in figure 16. Examination of these static co- 
efficients indicates a degree of nonlinearity in botli the 
normal-force and pitching-moment curves which may be due 
to the possible causes mentioned previously. 
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(a) Magnitude. 

(b) Phase angle. 

Variation of normal-force coefficient with reduced frequency for two-dimensional 65 AO 10 airfoil oscillating in pitch i 

the midchord at Mach number of 0.35. 
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(a) Magnitude. 

(b) Phase angle. 

Figure 10. — Variation of pitching-moment coefficient with reduced frequency for two-dimensional 65 AO 10 airfoil oscillating in 

pitch about the midchord at Mach number of 0.35. 


Oscillating normal-force coefficient. 
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(a) Magnitude. 

(b) Phase angle. 

Figure 11. — Variation of normal-force coefficient with reduced frequency for two-dimensional 65A010 airfoil oscillating in pitch about 

the midchord at Mach number of 0.7. 


Oscillating moment coefficient, 
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(b) Phase angle. 

Figure 12. — Variation of pitching-moment coefficient with reduced frequency for two-dimensional 65 AO 10 airfoil oscillating in 

pitch about the midchord at Mach number of 0.7. 


418788—57 3 
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Figure 13.— Variation of damping-moment coefficient with reduced frequency for two-dimensional 65A010 airfoil oscillating in pitch 

about the midchord at Mach number of 0.35. 
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Figure 14. — Variation of damping-moment coefficient with reduced frequency for two-dimensional 65A0I0 airfoil oscillating in pitch 

about the midchord at Mach number of 0.7. 
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Figure 15. — -Variation of static-normal-force and static-pitching- 
moment coefficients with angle of attack for two-dimensional 65 AO 10 
airfoil at Mach number of 0.35. 



Figure 16. — Variation of static-normal-force and static-pitching- 
moment coefficients with angle of attack for two-dimensional 65 AO 10 
airfoil at Mach number of 0.7. 

HIGH ANGLES OF ATTACK 

The normal-force and pitching-moment coefficients meas- 
ured at angles of attack higher than a=±° (figs. 9 to 12) 
differ appreciably from the coefficients measured at lower 
angles of attack. At zero frequency and high mean angles 
of attack, the magnitude of the change in normal force or 
pitching moment with change in angle of attack is consider- 
ably reduced from that obtained at small mean angles. In 


addition, in some cases the sign of the change in force or 
moment is reversed as indicated by values of the phase angle 
occurring at 180°. At frequencies other than zero, however, 
the magnitude of the forces and moments in some cases are 
much higher than those observed at low angles. 

Inasmuch as the coefficients of normal force and pitching 
moment presented in figures 9 to 12 correspond to slope of the 
lift curve or moment curve, it may appear somewhat sur- 
prising to observe much higher values for angles of attack 
near the stall than occurred at small angles of attack. This 
result lias been obtained before. (See refs. 14 and 15, for 
example.) Although the mechanism of this phenomena has 
not been established, a few speculative remarks may be made. 
If it is reasoned that the airfoil with separated flow may be 
replaced by a new airfoil having dimensional characteristics 
determined by the boundary layer, it does not appear in- 
conceivable that this imaginary airfoil could, under some 
oscillating conditions, undergo large changes in camber 
during the oscillation which, in turn, could produce much 
larger changes in normal force and pitching moment than 
would be expected of a nondeforming airfoil. Jones (ref. 4) 
has incorporated this feature of a variable-camber airfoil 
with some success in an attempt to calculate oscillating air 
forces at high angles of attack. However, the establishment 
of the time behavior of this imaginary airfoil would be a 
complicated problem. 

In addition to the large effects of increased angle of attack 
on the magnitudes of normal force and pitching moment, 
there were also large effects on the phase angles. The effect 
of angle of attack on the moment phase angles is such that 
under some conditions the aerodynamic-damping moment is 
unstable. This result is illustrated in figures 13 and 14 where 
the imaginary or damping component of the pitching- 
moment coefficient has been plotted as a function of reduced 
frequency. The existence of regions of unstable damping 
implies the possibility of single-degree-of-freedom flutter in 
the torsion mode. 

Examination of figures 13 and 14 indicates that the effect 
of increasing the Mach number from 0.35 to 0.7 was primarily 
to decrease the initial mean angle of attack at which un- 
stable damping moments occurred. This result, presum- 
ably, was due to larger adverse pressure gradients which 
caused boundary-layer effects. 

The existence of the large regions of angle of attack and 
reduced frequency at which unstable damping in pitch 
occurs is in direct contrast to the measurements of damping 
in the bending mode of finite-span wings which will be dis- 
cussed subsequently. 

FINITE-SPAN WINGS OSCILLATING IN FIRST 
BENDING MODE 

TECHNIQUE 

The aerodynamic damping for finite-span wings oscillating 
in the first bending mode was measured by use of a tech- 
nique based on measurements of the power required to main- 
tain oscillation at the natural frequency of the system. In 
order to eliminate from the measurements the presence of 
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relatively large random forces which acted on the wings at 
high angles of attack, a wattmeter was used as an indicator of 
the power required to maintain oscillation. It may be re- 
called that a wattmeter is a phase -sensitive device; and, since 
the phase of the random forces relative to the position of the 
wing is just as likely to cause a momentary increase in the 
power reading as it is to cause a momentary decrease, there 
will be no net effect of the random forces on the average 
reading of the wattmeter. A more detailed description of 
the technique follows. 

DESCRIPTION OF APPARATUS 

Wind tunnel. — The Langley 2- by 4-foot flutter research 
tunnel, described previously, was also used for the tests of 
the finite-span models. 

Models and oscillator. — The two semispan models used for 
the bending measurements differed primarily in airfoil thick- 
ness. The 10-percent-thick wing was rectangular in plan 
form and had NACA 65A010 airfoil sections. The wing had 
a chord of 8 inches, a semispan of 20 inches, and thus an 
aspect ratio of 5.0. The strength and stiffness of the 10- 
percent wing were concentrated in an aluminum-alloy core 
/32 inch thick, the airfoil being filled out by balsa glued to the 
core with the grain oriented vertically to the plane of the 
wing. The 3-percent-thick wing was tapered in plan form 
and had NACA 16-003 airfoil sections. This wing had a 
root chord of 10 inches, a tip chord of 5 inches, and thus a 
taper ratio of 0.5. The semispan of 16 inches resulted in an 
aspect ratio of 4.27. The 3-percent wing was constructed of 
solid aluminum alloy. Both the 10-percent wing and the 
3-percent wing had zero sweep of the midchord axis. 

The finite-span models were clamped at the root and os- 
cillated in, essentially, the first bending mode by a wire 
attached to the wing and to an electromagnetic shaker as 
shown schematically in figure 17. Since it was necessary 
to maintain a net tension in the wire between the shaker and 
the wing, in order for the shaker to apply a sinusoidal force 
to the wing, an external cantilever spring was installed 
between the shaker and the wing. The net effect of the 
additional stiffness and mass of the external wire-spring- 
shaker-coil system was to increase the first natural frequency 
of the 10-percent wing from 11.3 cps to 26 cps. The in- 
crease for the 3-percent wing was from 43 cps to 50 cps. 
The effect of the additional stiffness and mass on the oscilla- 
tory deflection modes for the two wings is illustrated in 
figure 18, where the shapes of the fundamental bending 
mode calculated for the wings as simple cantilevers are 
compared with the measured mode shapes of the wings 
oscillating at the first natural frequency of the combined 
system. The deflection modes were measured by use of a 
time-exposure photographic technique. (See ref. 16.) 

Instrumentation. — Both the 10-percent wing and the 3- 
percent wing were equipped with resistance-wire-strain-gage 
bridges arranged to indicate strains in bending and torsion. 
The bending gages were used as a qualitative measure of 
amplitude of oscillation. 

A block diagram of the instrumentation used in conjunc- 
tion with the electromagnetic shaker is shown in figure 19. 



Figure 17. — Schematic drawing of test setup for three-dimensional 

measurements. 

A variable-frequency voltage was supplied to the shaker by 
a low-frequency oscillator and a power amplifier. The 
voltage across the coil of the shaker was measured by use of 
a vacuum-tube voltmeter. The current through the coil 
was measured by use of a vacuum-tube voltmeter shunted 
with a 1-ohm precision resistor. The power consumed by 
the shaker was measured by use of a low-power-factor 
wattmeter. 

A photograph of the research instrument console used in 
the Langley 2- by 4-foot flutter research tunnel is shown in 
figure 20. 
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Figure 18. — Comparison of actual bending modes for three- 
dimensional tests with calculated fundamental bending modes for 
same wings as ideal cantileveis. 



Figure 19. — Block diagram of instrumentation used for three- 
dimensional tests. 



Figure 20. — Photograph of research instrument console used in 
Langley 2- by 4-foot flutter research tunnel. 
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PROCEDURE 

The procedures used in determining the aerodynamic 
damping in the bonding mode were essentially the same for 
the two finite-span wings except that the data for the 10- 
percent wing were obtained at fixed velocities with angle of 
attack as the primary variable, whereas the data for the 
3-percent wing were obtained at constant angles of attack 
with velocity as the primary variable. 

The procedure used in the tests of the 10-percent wing 
was as follows: After the desired velocity had been set, the 
mean angle of attack was adjusted by a motor-driven angle- 
of-attack changing and clamping device. The frequency 
of the oscillator was varied until the system was tuned to 
the frequency of maximum response. The pov r er to the 
electromagnetic shaker was adjusted to bring the amplitude 
of oscillation to the desired level as indicated by a mean- 
square root-bending-strain indicator. The desired voltage, 
current, and power were then observed. 

The preceding procedure was repeated in increments of 
mean angle of attack of about 2° up to 20°. Because of the 
occurrence of stall flutter in a torsion mode, the limiting 
angle of attack at the higher velocities was less than 20°. 

Essentially, the same procedure was followed for the 3- 
percent wing except that, at each angle of attack, measure- 



(a) 3-percent-thick wing (b) 10-percent-thick wing. 

Figure 21.- — -Variation of Mach number, Reynolds number, and tip angularity with reduced velocity for finite-span models. 


ments were made at several velocities before the angle was 
changed . 

The aerodynamic-damping coefficients for bending os- 
cillations of the 10-percent and the 3-percent wings were 
determined at several values of reduced frequency by varying 
the fiow r velocity. The tests were conducted at atmospheric 
pressure in air. The variations of Mach number, Reynolds 
number, and tip angularity with reduced velocity (the 
reciprocal of reduced frequency) for the two wings are shown 
in figure 21. 

DATA REDUCTION 

For clarity, a somewhat simplified presentation of the 
relationships involved in reducing the data is developed 
here. A more detailed presentation is made in appendix D. 

The actual distributed system used for the three-dimen- 
sional measurements is replaced by an equivalent single- 
degree-of-freedom system which has the following equation 
of motion: 

M e h+(C a +C 8 )h+KJi=F(t) (1) 

where h is the bending deflection of the wing and the coef- 
ficients M e and K c represent the generalized mass and stiff- 
ness, respectively. The coefficients C a and C s refer to the 
aerodynamic and structural damping, respectively. The 
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external force F(t) is again composed of three parts: a 
steady-state force associated with the mean angle of attack 
a, a sinusoidal force F 0 , and a random force F R . The 
steady-state force leads only to a static deflection of the 
beam and consequently will not be treated further. 

Equation (1) may be rewritten as 

M e (h 0 -{-h R ) + (C a -\-C s ) (h 0 -\rh R ) -\-K e (h 0 -{-h R )=F Q -{-I i R ( 2 ) 

where the subscripts 0 and R refer to quantities related to 
the forces F 0 and F R , respectively. 

Multiplying equation (2) by F 0 and taking time averages 
yields 

Kj^ 0 =T?+f*Ft 

where bars over the quantities refer to a time average over 
an arbitrary period of time T. However, 

MJirF o= 0 
( C a + Cs ) h R Fq = 0 
K e h R F 0 = 0 

FrF 0 =0 


each test so that it could be subtracted from the values 
measured during the test. 

Of course, certain corrections are required, and certain 
approximations are involved in converting the measured 
electrical quantities to the desired damping coefficient. For 
example, the electrical -power losses in the shaker coil must 
be extracted as a tare. In addition, the random force F R 
produces random velocities of the coil h R which, in turn, in- 
duce correlated voltages and currents in the electrical system. 
These currents are small, however, if the driving coil is 
excited by a high impedance generator. 

The ratios of aerodynamic damping to critical damping 
obtained in this manner apply only to the particular me- 
chanical systems used in these experiments. In order to 
make the damping measurements of more general value, the 
ratios of aerodynamic damping to critical damping were 
used to determine an effective section derivative F e having 
the same form as Theodorsen’s circulation function F(k). 
(See ref. 17.) The relationship between the damping ratio 
and the effective aerodynamic-damping coefficient is de- 
veloped in appendix D. The derivative F e is defined as 


as T — ^ go because the time average of the product of two 
uncorrelated quantities is zero. Also, at the natural fre- 
quency of the system, the sum of kinetic and potential 
forces must be zero; this fact leads to the condition that 

MJiqFq + K e h 0 Fo = 0 

Therefore, 

(Ca+C s )J^o=FV 

or 



KF, o 


Critical damping of the system is defined as 

C cr =2^M e 

where is the first natural frequency of the system and is 
related to the stiffness and mass by the following equation: 


Then, the ratio of aerodynamic damping to critical damp- 
ing may be defined as 

^± = IL 1 

( ' Cr h 0 F 0 Ccr 

The aerodynamic-damping ratio can be determined by 
measuring the current through the shaker coil which is pro- 
portional to F 0 and by measuring the power consumed by 
the system which is proportional to the product hF 0 . The 
structural-damping ratio was measured in still air before 


where k is the reduced frequency based on the chord at the 
three-quarter semispan and K e is a mass-ratio parameter 
shown in appendix I). Values of F e determined in this 
manner represent an effective or weighted section derivative 
which may be used in a strip-analysis estimation of damp- 
ing in the bending mode for wings having geometric prop- 
erties similar to the two tested in this investigation. 

RESULTS AND DISCUSSION 

The measured ratios of aerodynamic damping to critical 
damping in the bending mode for the 10-percent wing are 
shown in figure 22 as functions of angle of attack for various 
values of reduced frequency. The damping ratios for the 
3-percent wing are shown in figure 23 as functions of reduced 
velocity for various angles of attack. The effective aero- 
dynamic-damping coefficients in the bending mode F e for 
both the 10-percent and the 3-percent wings are shown in 
figures 24 and 25 as functions of reduced frequency at various 
angles of attack. 

LOW ANGLES OF ATTACK 

As a matter of reference, Theodorsen’s theoretical two- 
dimensional circulation function F(k) (ref. 17) is shown in 
figures 24 and 25 which are plots of the measured effective 
aerodynamic-damping coefficients for bending oscillations as 
functions of reduced frequency. As might be expected, the 
data for these relatively low-aspect-ratio wings at the lower 
mean angles of attack generally fall considerably below the 
values of F e indicated by the two-dimensional theory. 

The data for the 10-percent and the 3-percent wings are 
generally similar at the 1ow t angles of attack, except at low* 
values of reduced frequency. 


Co/C, 
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Figure 22. — Variation of aerodynamic-damping ratio for bending oscillations with angle of attack for 10-percent-thick wing. 
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Figure 23. — Variation of aerodynamic-damping ratio for bending oscillations with reduced velocity for 3-percent-thick wing. 


Damping coefficient, 
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Figure 24. — Effective aerodynamic-damping coefficient for bending oscillations for 10-percent-thiek wing. 


t* 
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Reduced frequency, k 


Figure 25. — Effective aerodynamic-damping coefficient for bending oscillations for 3-percent-thick wing. 
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HIGH ANGLES OF ATTACK 

In general, the outstanding characteristic of the damping 
measurements for bending oscillations presented in figures 
22 to 25 is the increased positive damping which occurred at 
the higher mean angles of attack. Within the range of angle of 
attack and reduced frequency covered by these experiments, 
the damping in the bending mode remained positive for all 
conditions; and at some conditions, the damping at high 
mean angles of attack was as much as about three times the 
damping at zero angle of attack. 

This increase in damping of the bending mode at high 
mean angles of attack may be related to the increase in the 
magnitude of normal-force and pitching-moment coefficients 
which was observed for the two-dimensional wing oscillating 
in pitch. However, for the pitching motion the increase in 
magnitudes was accompanied by changes in phase angle 
which produced negative damping moments. This difference 
in damping characteristics is not inconsistent with the results 
obtained for the bending measurements since the phase 
relationships which determine aerodynamic damping are 
very different for the two types of motion. Generally 
speaking, in the frequency range of practical interest, the 
predominant source of aerodynamic force or moment for a 
pitching wing is associated with the position terms, whereas in 
the case of bending oscillations the predominant source is 
associated with the velocity terms. 

CONCLUSIONS 

The magnitude and phase of the fundamental components 
of normal force and pitching moment acting on an NACA 


65A010 airfoil oscillating in pitch about the midchord have 
been measured at both high and low mean angles of attack 
at Mach numbers of 0.35 and 0.7. In addition, the aerody- 
namic damping in, essentially, the first bending mode has 
been measured for two finite-span wings over a range of mean 
angles of attack and reduced frequency. Examination of the 
results of this investigation indicate the following significant 
conclusions : 

1. The magnitudes of oscillatory normal-force and pitch- 
ing-moment coefficients for a two-dimensional wing oscillat- 
ing in pitch about the midchord at high mean angles of attack 
are much larger at some conditions than at low mean angles 
of attack. 

2. Large regions of angle of attack and reduced frequency 
exist wherein a single-degree-of-freedom torsion flutter is 
possible because of unstable aerodynamic-damping moments 
for pitching oscillations. 

3. The effect of increasing Mach number from 0.35 to 0.7 
was to decrease the initial angle of attack at which unstable 
damping occurred. 

4. The aerodynamic damping for bending oscillations was 
generally higher at high angles of attack than at low angles 
of attack, and no unstable damping was encountered over 
the range of configurations, angle of attack, and reduced 
frequency covered by these experiments. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 21, 1956. 


APPENDIX A 

effects of leading-edge conditions on aerodynamic damping for airfoil oscillating in pitch 


The effects of two types of disturbances placed at the 
leading edge of an airfoil oscillating in pitch were determined 
at a Mach number of 0.35 and angles of attack of 8° and 12°. 
The two types of disturbances which were chosen for study 
were a leading-edge roughness which consisted of a sprinkling 
of No. 120 carborundum grains over the leading 4 percent of 
the chord and a tripper-wire configuration which consisted 
of a wire about 0.4 percent of the chord in diameter located 
on the upper surface at 5 percent of the chord rearward of 
the leading edge. The magnitude and phase angle of the 
fundamental component of pitching moment were measured. 
The damping-moment coefficients obtained at the two mean 
angles of attack tested are shown in figure 26 along with the 
values obtained for the clean airfoil. 

At a mean angle of attack of 8°, the effect of the changes in 
leading-edge condition is pronounced. The damping co- 
efficients for the leading-edge-roughness tests reach con- 
siderably larger unstable values and the crossover point 
between stable and unstable values occurs at a much higher 


value of reduced frequency. The data for the tripper-wire 
tests at a — S°, on the other hand, indicate some improve- 
ment over the clean airfoil. The conventional static-pitch- 
ing-moment coefficients have been measured for the three 
airfoil conditions and are shown in figure 27. It appears that 
the addition of the tripper wire to the upper surface of the 
airfoil increased the angle of attack at which the maximum 
static-pitching-moment coefficient occurred from about 9° to 
about 10°, whereas the addition of the leading-edge rough- 
ness decreased the angle of maximum moment to about 8°. 
It appears, therefore, that the severity of unstable damping 
in this case depends upon the nearness of the mean angle of 
attack to the angle of attack at which maximum static 
pitching moment occurs. The damping moments measured 
at a = 12°, shown in figure 26, are relatively unaffected by the 
condition of the leading edge. This insensitivity, presum- 
ably, is due to the fact that at 5 = 12° the airfoil is well 
beyond the angle of maximum moment (as indicated in fig. 
27) for all three conditions. 
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(a) a = 8°. 

(b) a— 12°. 

Figure 26. — Variation of damping-moment coefficients with reduced frequency for 65A010 airfoil having three different leading-edge conditions. 


Moment coefficient, c m 
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Figure 27. — Variation of static-pitching-moment coefficient with angle of attack for 65 AO 10 airfoil having different leading-edge conditions. 
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APPENDIX B 

SOME OBSERVATIONS ON TIME HISTORIES OF INDIVIDUAL PRESSURES ON AN OSCILLATING AIRFOIL NEAR STALL 


Numerous investigators have attempted to formulate an 
analytical approach to the prediction of unsteady air forces 
for wings oscillating near the stall. Since this problem 
inherently involves the complexities of unsteady boundary- 
layer behavior, almost all the approaches have been em- 
pirical in the sense that experimentally observed phenomena 
have been utilized in formulating the necessary assumptions. 
So far, these analytical calculations have yielded satisfactory 
results for only limited ranges of the various parameters. 
One contributing cause of this lack of a satisfactory analyti- 
cal approach is that the phenomena involved in the process 
have not been completely described experimentally. In 
order to contribute further to the store of knowledge con- 
cerning oscillating air forces for stalled wings, some addi- 
tional information obtained during the course of the investi- 
gation of the NACA 65A010 airfoil pitching about the 
midchord is included. 

At a Mach number of 0.35 and Reynolds number of 
5.3X10 6 , a unique condition was encountered at angles of 
attack near 8°. A series of short samples of oscillographic 
records of the individual pressures is reproduced in figure 28. 
At an angle of attack of 6.8° and at zero frequency, the 
incremental pressures over the wing were relatively undis- 
turbed except for some small variations which occurred near 


the trailing edge. As the angle of attack was slowly in- 
creased, intermittent pressure fluctuations occurred near 
the leading edge at angles of attack above about 7.6°. This 
condition is illustrated in figure 28 (b) for an angle of attack 
of 8° and zero frequency. As the angle of attack was further 
increased, all the pressure gages indicated a continuous 
random fluctuation which is illustrated in figure 28 (c) for 
an angle of attack of 9.2°. 

When the airfoil was oscillated in pitch about the mid- 
chord at a mean angle of attack of 8°, this pattern of flow 
breakdown indicated by the random fluctuations in pressure 
was maintained in the sense that the flow periodically broke 
down as the angle of attack was periodically increased. The 
primary effect of the oscillation appeared to be a delay in 
the instantaneous angle of attack at which the flow break- 
down began. At a reduced frequency of 0.073, the gage 
nearest the leading edge indicated a disturbance commencing 
at an angle of attack of about 7.8° which represents a time 
delay of about 6 milliseconds beyond the time at which the 
airfoil reached 7.6° (the angle at which flow breakdown 
began for static conditions). At a reduced frequency of 
0.270, the angle of attack at which flow fluctuation began 
was delayed to about 8.5° which corresponds to a time delay 
of about 8 milliseconds. 
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Figure 28. — Time-history records of individual pressures of 65A010 airfoil oscillating in pitch at Mach number of 0.35. Positive pressure and 

positive angle of attack is up. 
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MEASUREMENT OF OSCILLATING AERODYNAMIC FORCES FOR VARIOUS MEAN ANGLES OF ATTACK 


31 


APPENDIX C 

SOME REMARKS CONCERNING NONLINEARITIES IN AERODYNAMIC FORCES 


The possibility of nonlinear aerodynamic effects in stall 
flutter and buffeting cannot be ignored. For linear systems, 
once the instability point has been passed, a divergent oscil- 
lation occurs. In stall flutter, however, it has been observed 
that the oscillations will usually reach an equilibrium ampli- 
tude for fixed conditions. Sisto (ref. 1) has pointed out 
that any attempt to estimate amplitudes of stall flutter 
must consider the nonlinear aspects of the phenomena. In 
buffeting studies of the type proposed by Liepmann (refs. 
5 and 6), it is necessary to assume that the system is linear. 
Thus, the problem of determining the nonlinear character- 
istics of oscillating air forces for stalled wings is important. 

In the present investigation, a few limited studies of the 
nonlinearity of the oscillatory air forces have been made. 
Two different techniques were employed. For the 10- 
percent-thick wing oscillating in the first bending mode, the 
aerodynamic damping was measured as a function of ampli- 
tude of oscillation. For the two-dimensional wing, the 
oscillatory normal forces and pitching moments at several 
conditions were examined for harmonic content. 

If a wing oscillates in simple harmonic motion, nonlinear- 
ities in the aerodynamic forces should cause higher har- 
monics to appear in the oscillatory air forces. Since the 
amplitude of oscillation of the two-dimensional wing could 
not easily be varied, the harmonic-content method of study- 
ing the nonlinearities was adopted. At several conditions 
of angle of attack and frequency, the electrical signals pro- 
portional to normal force or pitching moment were scanned 
by tuning the narrow band-pass filter through its frequency 
range. No significant harmonic content was noted for any 
of the conditions examined. Figure 6 is representative of 
the results obtained. Since the amplitude of oscillation 
used in the two-dimensional tests was relatively small 
(nominally ±1.2°), it is possible that the nonlinearities 
which may have existed were not sufficiently large to be 
detected by the method used. 


The results obtained for the 10-percent-thick wing oscil- 
lating in the bending mode at various tip amplitudes are 
shown in figure 29 where the aerodynamic-damping ratio is 
plotted as a function of tip amplitude for various angles of 
attack at a reduced frequency of about 0.35. The range of 
tip amplitudes covered in this experiment represents a tip 
angularity from ±0.66° to ±2.3°. The data do not indi- 
cate any particular trend of aerodynamic damping with 
amplitude; however, there are differences of about 25 per- 
cent at various amplitudes. The significance of differences 
of this size would, of course, depend upon the nature of the 
calculation in which the damping values might be used. 



Figure 29. — Variation of aerodynamic-damping ratio for 10-percent- 
thick wing with amplitude of oscillation. fc«0.35. 
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APPENDIX D 


DEVELOPMENT OF DATA- REDUCTION RELATIONS FOR FINITE-SPAN WINGS OSCILLATING IN FIRST BENDING MODE 


The more important features of the method employed to 
deduce the desired nondimensional coefficients of aero- 
dynamic damping for the finite-span wings oscillating in the 
first bending mode were presented in the text. The purpose 
of this appendix is to present the development of some of the 
relations used in reducing the data in somewhat more detail. 
Some of the features of the method which are presented in 
the text (such as the manner of accounting for the random 
force) are not included. 

The system is a cantilever beam which has attached to it 
at the spanwise coordinate y c an external mass M c and 
external spring K c . The external mass and spring act on 
the beam through the mean angle of attack a. as indicated in 
sketches 1 and 2. The S3 r stem is excited by a vertical 




sinusoidal force F 0 e iu}t which also acts upon the beam through 
the angle a. It is assumed that the motion of the beam is 
confined to one degree of freedom, namely, the first bending 
mode of the S3^stem which is defined as h=a 1 h l (y) where h is 
the bending deflection of the beam normal to the chord, a x is 
the tip deflection, and h x (y) is the first-bending-mode shape 
expressed in terms of a unit tip deflection. The following 
development is similar to that in reference 18. 

The potential energ3^ of this system is 

E P =~j o EI ^^2 Mi (*/)] dy+\ K c cos aj (Dl) 

The effective stiffness of the beam in the first bending mode 
is defined as 

K h =£EI^My))ydy 


The kinetic energy of the system is 

i 2 j o ^[^i(y)] 2 dy+^di 2 M c [hi(y c ) cos a] 2 (D2) 

The effective mass of the beam in the first bending mode is 
defined as 

Mhl= j 0 m[h x (y)] 2 dy 
and equation (D2) becomes 

E t=\ dx 2 M c [h x (y c ) cos a] 2 

=1 d?M e 

If the assumption of viscous damping for the system is made, 
the structural damping is described by the dissipation func- 
tion which is defined as 

D =\ dS^CMyWdy 

or 

D=\iSC s 

The virtual work 8 W done by the external forces acting,on 
the S3"stem when the system is displaced through a virtual 
displacement 8a x is 


8\V= — TrpVF e cos aj 8aicdih l 2 (y)dy-\-F 0 8aihi(y c ) 


cos a 


The generalized force is then 


=—irpVdiF e cos clhxiy^dy+FJiiiyc) cos a 


The aerodynamic forces associated with acceleration and 
displacement have been ignored in this development. The 
experimental technique, however, accounts for these forces 
in that the measured natural frequency for each condition 
was used in reducing the data. B3" use of Lagrange’s 
equation of motion (ref. 18) it can be shown that 


Thus, equation (Dl) becomes 

E P =\ 2 K hx -\-^ ai 2 K c [hi(y c ) cos S] 2 


where 


M e di~\- (C s -\-Ca)di-\-K e (ii — F 0 hi(y c ) cos a 
C a =TpVF e cos c[h x (y)] 2 dy 



It is desirable to determine the ratio of aerodynamic 
damping to critical damping of the system. Critical 
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damping is defined as 


Since 

then 

and 


C cr = 2yKM 

M e o>i 1 2 =K e 

C cr =2u 1 M e 

c TT pVF e cos aj^ c[hi(y)] 2 4 dy 
C c r 2a>! M c 


This expression may be rearranged to solve for a more 
general nondimensional aerodynamic-damping coefficient F c 


F,: 


tt 2u l M e 

^ cr 


7 rpv COS a j. c[hi(y)]' 2 dy 

and in terms of conventional flut ter parameters 


F f 


^ cr 

K e COS a 


where 


and 



C 

In order to solve for either - - or F e , the natural frequency 

C cr 

oi\ and the generalized mass M e must be known. 

As mentioned previously, coj was determined experi- 
mentally at each test condition — a procedure which auto- 
matically accounts for aerod 3 r namic-st iff ness effects. The 
generalized mass of the system was also determined experi- 
mentally by measuring the change in natural frequency 
produced by attaching a known concentrated mass to the 
shaker coil. Tiie natural frequency of the system without 
the additional mass is 



and the natural frequency of the system with the added 
mass is 

2 _ K e 

M,+AM c My e ) cos a] 2 

where A M c is the added mass. From these expressions for 
the natural frequency, it follows that 


M P = 


AM c [hi(y c ) cos a] 


«ia S 6 7 8 9 


*1 


Since the determinations of M e were made in still air rather 
than in a vacuum, the values of generalized mass obtained 
would correctly include the apparent mass of the air sur- 
rounding the wing. 
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